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Abstract. In this paper we study the extended atmosphere of the late-type supergiant a Orionis. Infrared spec- 
troscopy of red supergiants reveals strong molecular bands, some of which do not originate in the photosphere 
but in a cooler layer of molecular material above it. Lately, these layers have been spatially resolved by near 
and mid-IR interferometry. In this paper, we try to reconcile the IR interferometric and ISO-SWS spectroscopic 
results on a Orionis with a thorough modelling of the photosphere, molecular layer(s) and dust shell. From the 
ISO and near-IR interferometric observations, we find that a Orionis has only a very low density water layer 
close above the photosphere. However, mid-IR interferometric observations and a narrow-slit N-band spectrum 
suggest much larger extra-photospheric opacity close to the photosphere at those wavelengths, even when taking 
into account the detached dust shell. We argue that this cannot be due to the water layer, and that another 
source of mid-IR opacity must be present. We show that this opacity source is probably neither molecular nor 
chromospheric. Rather, we present amorphous alumina (AI2O3) as the best candidate and discuss this hypothesis 
in the framework of dust-condensation scenarios. 



Key words. Techniques: high angular resolution - Techniques: spectroscopic - Stars: individual: a Orionis - Stars: 
atmospheres - Stars: supergiants - Stars: circumstellar matter 



1. Introduction 

\ It is now becoming clear that late- type supergiant stars, 
like their less massive AGB counterparts, are also embed- 
ded in a circumstellar environment ( CSE) of molecular 
layer s gaseous outflow and dust (e.g. iRichards k. Yates] 
Il998t iTsuiil l2000ri IPerrin et alJl2005|) . What sets them 
apart, however, are the low amplitude pulsations, rel- 
atively high effective temperatures and sometimes the 
presence of a chromosphere. It is therefore questionable 
whether the mechanism driving their mass-loss is similar 
to that of AGB stars, i.e. a complex interplay between 
pulsations, shock waves, molecu l ar opacity and dus t con- 
densation l|Fleischer et al.lll995l iHofner et al"]l2003l and 
references therein). 

a Orionis and fi Cephei play an important role in 
this area of research, since it was in these stars that 
molecular layers (containing water) around late-type su- 
pergiants were detected for the first time: ITsuiil |l978|) 
found u Cephei to be 0.5 mag brighter between 5 and 



8 /im due to emis sion by hot O in the circumstellar en- 
r. lTsuiill|2000a|i 



vironmcnt. Later, 



(2000a) also found emission lines 
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of water in the ISO-SWS spectrum of fi Cephei and deter- 
mined temperature, location and column density of this 
so-called molsphere. In a r e-analy sis of old Stratoscope 
II data of a Orionis, ITsuiil l)2000bj) identified unexpected 
absorption lines as due to non-photospheric water. 

Their large apparent size makes both supergiants good 
candidates for optical/IR interferometry. In recent years, 
the presence of these molecular layers was confirmed with 
optical/IR interferometry, allowing the determination of 
their size, temp erature and optica l dept h at the wave- 
lengths studied l|Perrin et alJl2004al 12005]) . 

Generally, modelling attempts of the molecular ex- 
tended atmospheres of these supergiants have targeted ei- 
ther spectroscopic data or interferometric data, but not 
both simult aneously. The o ne exception is a study of 
a Orionis by lOhnakal ll2004blb in which mid-IR interfero- 
metric and spectroscopic data are modelled using a black- 
body+molecular layer approximation. 

In this paper we present a model for a Orionis, which 
fits near- to mid-IR interferometric and spectroscopic 
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data simultaneously. Other crucial difference with pre- 
vious works is that we take the photospheric molecular 
bands into account and cover a larger wavelength regime. 
In section |21 we present the IR data collected from the lit- 
erature and the supplementary data used in our analysis. 
Sect. describes the models for photosphere, molecular 
layers and dust shell used to interpret the data. In Sect. 01 
we compare the models with data for a Orionis and dis- 
cuss the results. Finally, in Sect. El we conclude and look 
ahead. 



2. The observations 

2.1. Spectroscopy 

a Orionis was observed with the ISO-SWS (Infrared Space 
Observatory, 1995-1998, Short Wavelength Spectrometer, 
2.38 — 45.2 //m at a spectral resolving power of up to 
R - 3800) on MJD= 50729 with the AOT01 template 
at speed 4. AOT01 refers to a single full-wavelength up- 
down scan for each aper ture with four poss ible scan speeds 
at degraded resolution ( Leech et al1l2002() . Speed 4 is the 
slowest scanning velocity, yielding a spectral resolution of 
800 - 1500. 

In total, the SWS uses 4 detector arrays associ- 
ated with the SW (Short Wavelength) and LW (Long 
Wavelength) gratings (SW: 2.38- 12.0 /mi and LW: 12.0- 
45.2 /im), with 12 elements each. Every 12-element array 
observes one wavelength-band, with band 1 ranging from 
2.38 - 4.08 /im, band 2 from 4.08 - 12.0 /im, band 3 from 
12.0 - 29.0 /im and band 4 from 29.0 - 45.2 /im (a def- 
inition of all sub bands is given in Table ^) . The flux is 
measured by integrating the photocurrent produced in the 
photoconductors on an integrating capacitance on which 
the light of a certain wavelength is collected. The volt- 
age over the capacitance is read out non-destructively at 
24Hz during an integration interval. An integration inter- 
val typically lasts for one or two seconds. At the end of 
the integration interval the capacitor is discharged to start 
a new integration. The slope of the non-destructive read- 
outs of the voltage over the integrating capacitance is a 
measure of the flux falling onto the detector. 

The data reduction procedure we follow uses the IA 
(Interactive Analysis) tools also used in the latest ver- 
sion of the pipeline (OLP10.1, OLP = Off-line Processing), 
complemented with a manual removal of glitches, bad de- 
tector data points and scan jumps. For a detai l ed de- 
scri ption of these procedures, we refer to ICamil ll2002t) 
and I Van Malderen et alJ l)2004|b Furthermore, we ch eck 
the correction for memory effects llKester et al.l l20ol by 
comparing both up and down scans with the final spec- 
trum. Joining of the sub bands was done by means of the 
overlap between the different bands, and starting from 
band Id which is believed to have the best absolute flux 
calibration. For a detailed description of the need for this 
procedure (problematic dark current subtrac t ion, p ointing 
errors, etc), we refer to lVan Malderen et al.l l)2004[) . 



Below, we discuss the peculiarities of the data reduc- 



tion. 



Presaturation a Orionis brightness (F„ ~ 20000 Jy at 
2.5 /im) proves to be especially challenging, mainly be- 
cause the dynamic range of the detectors was optimised 
for sources of a few to a few thousand janksy. While there 
is no sign of actual saturation of the detectors, i.e. there 
is no clear cut-off value, the integration ramp of some of 
the detectors shows a strong non-linear behaviour when 
the voltage across the integrating capacity reaches very 
high values. This flattening of the integration-ramp to- 
ward the end of the integration time causes an underesti- 
mation of the slope and therefore of the flux at that wave- 
length. A similar phenomenon has been seen by the SWS- 
Instrument-Dedicated-Team (SIDT) in a post-Helium ob- 
servation (after the satellite had run out of coolant), and 
is called "presaturation" . In the pipeline-reduced data, 
this presaturation causes an "absorption feature" at 3 /im 
which is not real. 

The solution to this problem of presaturation, pre- 
sented by the SIDT, is straightforward: only use that part 
of the integration-ramp which behaves linearly. In con- 
crete this comes down to using only the first half or equiv- 
alently the first 24 voltage measurements minus those re- 
jected for other reasons 1 . We first tested this procedure 
on the ISO-SWS primary calibrator a Bootis and found 
that it does not introduce any other artefacts. 

A comparison between the 2 spectra, the original one 
and a corrected one, after reduction, is shown in Fig.^ the 
spurious spectral feature around 3.0/im has disappeared. 



Memory effects Band 2 (Si:Ga detectors, 4.08 — 12 /im) 
is notoriously sensitive to memory effects, especially for 
sources with a flux above 1000 Jy such as a Orionis 
l)Kester et al.ll200ljb The detector 'remembers' a previ- 
ous illumination and therefore predicts too large a flux for 
the current observation. Since each sub band is scanned 
both from short to long wavelengths and vice versa, it 
is possible to estimate the strength of the memory ef- 
fects. Fig.[21shows up scan, down scan and final version of 
the a Orionis spectrum. Strong memory effects are indeed 
present. We must therefore keep in mind that the correc- 
tion may not be flawless from 4.08 to 4.3 /im and between 
4.8 and 5.5 /im. 



1 In the case of this AOT01 speed 4 observation, the inte- 
gration time is 2 seconds and the voltage is read out 24 times 
per second. Hence, we have an integration-ramp consisting of 
48 points. Although decreasing the number of data points on 
which the fit is based might sound like a loss in signal-to-noise, 
one should bear in mind that those points not used now, made 
a large contribution to the variance with respect to the fitted 
slope because of their non-linear behaviour. As a result, the 
error bars actually decrease significantly. 
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Fig. 1. The final band 1 spectra of a Orionis : either 
based on uncorrected data (black) or on data corrected for 
presaturation (grey) . The spurious spectral feature around 
3.0/im has disappeared. 
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Fig. 2. Up scan (green), down scan (red) and final spec- 
trum (black) in band 2a and 2b. Memory effects are strong 
at the blue side and at 5.4 /im. The final flux may not be 
estimated correctly from 4.08 to 4.3 /jm and between 4.8 
and 5.5 /im. We keep this in mind when comparing to the 
models. 

Absolute calibration Another peculiarity of the reduction 
of this observation concerns the multiplicative 2 factors we 
use to paste the different sub bands together. The general 
strategy is the following: band Id (3.02 — 3.52 /im) is be- 
lieved to have the best absolute flux calibration. Using the 
overlap between the different sub bands, it is then possi- 
ble to shift the other bands to the correct level, starting 
from those adjacent to Id and moving band- by-band to 
the red and blue. Usually the factors used are quite close 
to unity, with differences up to a few percent. However, 

2 for bright sources like a Orionis, the additive terms, such 
as the dark current, are negligible in comparison to the multi- 
plicative effects due to e.g. pointing errors, problems with the 
spectral response function, etc. 



for the spectrum of a Orionis, several factors differ signif- 
icantly from unity. We remark that scaling factors closer 
to unity would be possible if we assume the absolute cal- 
ibration of band Id to be too low (see Sect. I4.3fl . Also 
the memory effects at 4.1 /im could influence the scaling 
factor of band 2a by 5%. For completeness, we list these 
scaling factors in Tabled 

Table 1. Multiplicative factors used to scale the differ- 
ent sub bands of the a Orionis spectrum to the absolute 
flux level when the latter is either determined by the flux 
level in band Id (3rd column) or by requiring the average 
of the scaling factors to be equal to 1 (4th column). The 
large deviations from unity in the ld-referenced scaling 
suggest that the absolute flux in band Id might be under- 
estimated. Further evidence supporting this hypothesis is 
given in Sect. 14.31 The large jumps between some bands 
can be (partially) attributed to the use of three different 
filters and apertures and 4 different detector types. 



Band 


A [/im] 


Id 


mean=l 


la 


2.38 - 2.60 


0.97 


1.15 


lb 


2.60 - 3.02 


0.98 


1.16 


Id 


3.02 - 3.52 


1.00 


1.19 


le 


3.52 - 4.08 


0.97 


1.15 


2a 


4.08 - 5.30 


0.78 


0.93 


2b 


5.30 - 7.00 


0.80 


0.95 


2c 


7.00 - 12.0 


0.66 


0.79 


3a 


12.0 - 16.5 


0.71 


0.84 


3c 


16.5-19.5 


0.75 


0.89 


3d 


19.5 - 27.5 


0.76 


0.90 



2.2. Interferometry 

Literature data a Orionis has been studied extensively 
with interferometric techniques in the Near- and Mid-IR. 
In this work, we use 3 sets of IR observations that were al- 
ready published: (1) bro ad band FLUOR obs ervations at 
2.2 /im first presented in lPerrin et~aH ll2004al), (2) T ISIS 
L band observations at 3.8 urn jcii^giioii et al.l l2002) and 
(3) ISI narrow-ban d Mid-IR observations at 11.15 /im 
JWeiner et a,l .ll20T)l . 

The K band data cover molecular bands of H2O and 
CO, the L band data those of H2O and SiO. These 
molecules may very well be present in the molsphere 
around a Orionis. We note that, if the extra molec- 
ular layers are optically thin, both bandpasses remain 
photosphere-dominated and are therefore well suited for 
our analysis. 

The mid-IR visibilities not only sample the photo- 
sphere+molecular layers, but also the dust shell which 
emits strongly at 11 /im. In the case of a Orionis, the 
silicate emission of the dust shell originates in a region far 
from the central obj ect: the inner radi us, Ri n , is believed 
to be 0.5 arcsecond llSloan et al.lll993h or even 1 arcsec- 
ond on the sky IjDanchi et alJll994h . to be compared to 
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22 mas for the photospheric radius (|Perrin et alJl2004a|) . 
The silicate emission is therefore totally resolved at the 
spatial frequencies of interest for the study of the central 
object. If one knows the flux ratio / = p Btar , it is possible 
to renormalize the observed visibilities allowing an inde- 
pendent study of the central object without any other in- 
formation on the dust shell. For a Orionis, this flux ratio 
at 11. 15 urn is well-determined through other ISI observa- 
tions l|Pa,nchi et al]ll994l Ifiester et a,l.lll99ri ISudol et all 
Il999l) at very low spatial frequencies which sample the 
visibility curve at the point where the dust shell is being 
resolved. The derived values for the flux ratio / range from 
55 to 65 percent. 

New data a Orionis was observed in the N-band 
with the Mid-IR i nterferometric instr ument MIDI (e.g. 
[ Leinert et aT1l2003[) on the VLT-I (e.g. iGlindemann et all 
2000) in Science Demonstration Time (SDT) on November 
8 and 10, 2003. The interferometric data will be presented 
and analysed in a forthcoming paper, together with ad- 
ditional MIDI observations scheduled for autumn 2005 
which will cover a wider range in spatial frequencies. 



MIDI single-dish spectra As shown in IChesneau et alJ 
(|2005[) . it is also possible to obtain good quality N-band 
spectra from MIDI observations: for calibration purposes, 
single-dish spectra are observed with both UT telescopes 
separately just prior or after the interferometric obser- 
vations. We searched these spectrally dispersed 1-D im- 
ages (0.54 arcsec wide slit) for extended emission from 
the detached dust shell but found none. Spectra were con- 
structed by integrating the images in the spatial direc- 
tion with or without mask, which did not yield signif- 
icant changes in spectral shape. We use here the spec- 
tra constructed with mask because of their higher S/N. 
They were calibrated using the MIDI spectra of calibrator 
stars (HD 39400, K1.5IIb and 17 Mon, K4III) observed in 
between the Betelgeuse observations. The final averaged 
spectrum is shown in Fig. ED 

3. The modelling 

The aim of this investigation is to simultaneously model 
IR spectroscopic and interferometric data in order to get 
a better understanding of the physical processes at work 
between photosphere and dust shell. Since we expect both 
the extended atmosphere and the dust shell to be optically 
thin, it is crucial that we treat every source of radiation 
properly. We can for example not neglect the strong molec- 
ular bands already present in the photospheric spectrum 
of a late- type supergiant: clearly not all molecular features 
seen in the IR originate in "extra" molecular layers. 

Hydrodynamical models of oxygen-rich late-type stars 
are not yet in a stage where they can provide an accurate 
reproduction of observations, and the general consensus 
is that we are still missing some fundamental knowledge 
on the major processes, such as the dust condensation se- 



quence (e.g. lTielensl [l990). Moreover, these models are 
currently built only for low-mass (Mira and C) stars, and 
their applicability to supergiants is not clear. Therefore, 
we choose not to use such a self-consistent model, but 
instead construct a semi-empirical model as described be- 
low. 

We will divide our model into 3 parts: (1) the hydro- 
static photosphere, (2) a region of extra molecular layers 
and (3) the dust shell. The approximation of the extended 
atmosphere, i.e. the part containing the extra molecular 
layers, with a discrete instead of a continuous density, tem- 
perature and composition distribution is supported by the 
first generation of hydrodynamical models for AGB stars 
which predic ts these distributions to be strongly ra dially 
peaked (e.g. IWoitke et al.lll993 iHofner et ablEofii. 

For each part, we have used the most up-to-date mod- 
els available. They are described below. 



3.1. The photosphere 

To model the photosphere of this supergiant, we use 
the sosmarcs code, ve rs ion May 1998, as devel oped by 
iGustafsson et all (|l975|) : IPlez et all (|l992l Il993j) at the 
university of Upssala, Sweden. This code is specifically 
developed for the modelling of cool evolved stars, i.e. a lot 
of effort was put into the molecular opacities, and they al- 
low for the computation of radiative transfer in a spherical 
geometry. The line lists — rele vant for the IR p art of the 
spec trum — used are those of iGoldma n et alJ 1119981) for 
OH. ICoorvitchl lll994h for CO. iLangfioff fc Bauschlichel 
l|l993h for SiO Mid those of IPartridge & Schwenkel i|l997t) 
for H 2 0. This code solves the radiative transfer equation 
in a hydrostatic LTE environment with an ALI 3 method 
fe.g. lNordlundlll984 . Opacities are treated with the opac- 
ity sampling (OS) technique at 153910 wavelength points, 
which guarantees a good sampling of both molecular and 
continuum opacity sources. This OS grid offers a good 
compromise between computational speed and accuracy 
of the atmospheric structure of the model. In the inter- 
pretation of the similarities and discrepancies between ob- 
served and synthetic spectrum, we make use of the study 
of the influence of the stellar p aram eters on a synthetic 
spectrum performed bv iDecir] l|2000h . The initial values 
of the model parameters (which are T e ff, logg, mass, the 
abundances of C, N, and O, the microturbulent velocity, 
£t and the metallicity [ Fe/HD were de t ermin ed in great 
detail for a Orionis bv lLambert et all l)l984l) . They are 
listed in Table 

For a Orionis, masses derived from theoreti- 
cal Hertzsprung-Russell diagrams range fro m 15 M Q 
( Lamb et al.lll976]ICloutman fe Whitakerlll980|) to 30 M 
( Stothers fc Chin! 19791). Adopting a distance of 131 par- 
sec fHiDDarcos. lPerrvman et alJll997)) and a photospheric 
angular diameter of 43 mas yields a stellar radius of about 
645 Rq . This result confirms the value for the gravity 
found bv lLambert et all l)l984|) . derived from the fact that 



ALI = Approximate Lambda Iterator 
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Table 2. The ste llar / atmospheric pa rameters of 
a Orionis, according to lLambert et all l)l984j) . for a tem- 
perature of 3800 K (their own estimate) and for a tem- 
perature of 3600 K (as suggested by other temperature 
determinations). 



parameter 


value 


value 


T e ff 


3800 K 


3600 K 


log5 


0.0 


0.0 


Mass 


15 M Q 


15 M Q 


e(C) 


8.41 


8.29 


e(N) 


8.62 


8.37 


e(O) 


8.77 


8.52 


€t 


4 km/s 


4 km/s 


[Fe/H] 


0.00 


0.00 



neutral and ionized lines should yield the same abun- 
dances. We stress the quality of this determination of the 
surface gravity because photospheric molecular bands, and 
especially th ose of water , are highly sensitive to the sur- 
face gravity l|Decml i2000b 

The molecules responsible for absorption in the ISO- 
SWS wavelength region are primarily CO, H2O, OH 
and SiO. The contribution of these molecules to the to- 
tal absorption, when adopting the stellar parameters of 
lLambert et all l)l984|) . is shown in Fig. 

3.2. The molecular layers 

Several approaches are possible when modelling mol- 
spheres. They range from pla ne-parallel slabs placed 
in front of the star (e.g. lYamamura et alJ jl£)99t 
iMatsuura et~a"fl l2002t ICamil 120021 IVan Malderenl 120031) 
over i nfinitesimally thin spherical shells flPCTrinet alJ 
20041)) to extended sph erical layers with temperature 
and density distributions l|Mennesson et alJl2002HOhnakal 
l2004albh. 

Slab models are not suitable for our purposes: the 
molecular layers are expected to be fairly optically thin, in 
which case the slab-approximation is far too rough (in a 
first order because these models do not take into account 
that the visible column density is twice as high next to 
the stellar disk as it is in front of it, and to a second order 
because they miss sphericity effects which are important 
for large scale heights). Both the resulting spectra and the 
resulting visibilities are therefore not realistic for optically 
thin 4 molspheres. 

On the other hand, the shocked nature of these layers 
most probably implies that they are geometrically quite 
thin and that therefore temperature and density gradients 



4 This does not take away their merit in the modelling of 
the thick molecular layers surrounding Mira-like stars, where 
the column density and consequently also the opacity are so 
large that the backside of the shell is not visible and the spatial 
intensity distribution is rather UD-like. 
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Fig. 3. Main contributions to the total photospheric ab- 
sorption i n the atmosphere mod el with the stellar param- 
eters from lLambert et all l|l984h at 3600 K. CO and SiO 
are clearly dominant, while there are only minor traces 
of absor ption by water. The relevant l ine lists used are 
those o f iGoldman et alJ dl998h for OH, iGoorvitcbl l|l994h 



for CO, Langhoff & Bauschlicher ( 1993) for SiO and those 
of Ames l|Partridge fc Schwenkelll997l) for H 2 0. 



within the shell can be neglected. Moreover, the observa- 
tional data are at the moment too limited (in both spectral 
and spatial resolution) to constrain these supplementary 
parameters. 

Consequently, we opt for isothermal, spherical layers 
with a finite extent but no density distribution, which are 
characterized by their temperature, composition, column 
density for each molecule, inner radius and outer radius. 
Fig. 0] visualizes the structure of the models used. 

3.2.1. Opacities 

Opacities were calculated for H2O, SiO, CO, and OH at 
temperatures ranging from 500 to 2500 K, with 100 K 
increments. 

For H2O, w e used the NASA AMES line list and par- 
tition function I Partridge fc Schwenkelir997lb which is at 
present the most complete list, including more than 300 
million lines. However, for computational reasons, we in- 
cluded only lines with \oggf — x@ > — 9? with g the 
ground level statistical weight, / the oscillator strength, 
9 = 5040/3500 and \ the excitation energy in eV. The 
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Fig. 4. The structure of the models presented here: a pho- 
tosphcric model which includes limb darkening and the 
relevant molecular opacity, surrounded by layers of possi- 
bly mixed composition, each with its temperature, density, 
inner radius and outer radius. For the silicate feature mod- 
elling, the resulting spe ctrum is fed into the du st radiative 
transfer code modust l|Bouwman et al.ll200o() . 



the resulting spectrum and the visibilities can be calcu- 
lated. For more det ails on the calcula tion of the radiative 
transfer, we refer to IVerhoelstl ((20051) . 

An example of such an intensity profile for a single 
layer around a limb-darkened central star is shown in 
Fig. El From the theoretical wavelength-dependent inten- 




0.0 L 
0.0 



0.5 1.0 1.5 

mpact parame;er [R ] 



Fig. 5. An example of an intensity profile at 2.9 /im for 
our MARCS photosphere, surrounded by a single layer with 
Ri n = 1.35R*, R out = 1.85R*, a temperature of 2300 K 
and a H2O column density of 5 x 10 20 cm -2 . The dip in 
the IP at ~ 1 R* is an artefact due to the assumption that 
the central star is optically thick also in the limb. 



consequences of neglecting the weaker lines on the final 
opac ity are below 1 percent fo r a spectral resolution of 
300 l|Van Malderen et all h()Q4\ . Differences between the 
available line lists for H 2 (e.g. AMES, HITEMP, scan) are 
significant. The hitemp line list is aimed at temperatures 
of about 1000 K, and might not be valid for t he very high 
tempe ratures of the water around Betelgeuse. lDecin et all 
( 2OO3J found the best reproduction of observed water spec- 
tra in the ISO-SWS spectra of M giants with the AMES line 
list, for which reason we used it here as well. 

CO, OH and SiO pr esent far les s difficulti es than water. 



We us ed th e line list of Goorvitchl lll994l). iQoldman et alJ 
l|l998j) and iLanghoff fc Bauschlicherl (|l993h respectively? 
These are the mo st complete l ists to date for astrophys- 
ical applications l|Decinl 12000'). The polynomial expan- 
sion of the appropriate partition function was taken from 
ISauval fc Tatuml l)l984l) . 

As microturbulent v elocity we use 3 km s _1 as sug- 
gested for cool giants bv I Aringer et alJ J2002I) . 



3.2.2. Radiative transfer 

We compute emerging intensities for a grid of linearly 
spaced impact parameters p (256 in total), i.e, we con- 
struct a 1-D intensity profile I(X,p), going from the center 
of the disk to the edge. From this intensity profile, both 



sity profile, we calculate the spectrum by integration over 
the emitting surface and the monochromatic visib ility by 
a Hankel transform (iHanburv Brown et al.Nl974[) of the 
intensity profile. This monochromatic visibility is then 
convolved with the bandpass of the observations for the 
comparison model vs. observations. For wide-band data 
close to and beyond the first null, bandwidth effects be- 
come important: because of the wavelength dependence of 
the spatial frequency at which we observe, one wide band 
visibility measurement actually covers a range in spatial 
frequencies, but is assigned to the spatial frequency cor- 
responding to the effe ctive wavelength of the filter. We 
follow the approach of IPerrin et alJ l|2004a|) and take this 
effect into account by averaging properly weighted squared 
visibilities. 

3.3. The dust shell 

The dust shell is modelled using the proprietary spherical 
radiative transfer code MODUST 

l|Bouwman et all l200fi iBouwmanl l200l|) . Under the 
constraint of radiative equilibrium, this code solves 
the monochromatic radiative transfer equation from 
UV/optical to millim etre wavelengths using a Fea utrier 
type solution method l)Feautrierll963 lMihalaslll978|) . The 
code allows to have several different dust components of 
various grain sizes and shapes. 
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4. Comparison with observations 

When confronting the synthetic spectrum based on our 
MARCS model atmosphere with the ISO-SWS spectrum 
of a Orionis (Fig. 0), it is clear that our model makes 
a good reproduction of the global spectral shape. This 
shape is quite different from that of a blackbody: clearly, 
the major sources of opacity have a strong influence on 
the spectral appearance in the IR, which is well repro- 
duced by the MARCS model atmosphere. However, several 
differences can be observed: 

— between 2.38 and 3.5 /mi the model strongly overesti- 
mates the flux, 

— from 4.0 to 4.2 /xm the flux is underestimated, 

— around 5 /xm, there are extra absorption features in 
the fSO-SWS spectrum, 

— also from 6 /xm onward there is extra absorption, and 

— from 8.5 /xm on- wards, the observed spectrum is dom- 
inated by silicate emission. 

The latter property clearly divides the spectrum in a 
part without dust contribution, and a part which will re- 
quire a modelling of the dusty outflow. We are therefore 
motivated to separate at first the analysis of the near-lR 
part of the spectrum from that of the mid- to far-IR wave- 
length regime. 

4.1. The near IR 

We start with a study of the near-IR wavelength regime, 
where photospheric emission dominates, but several dis- 
crepancies require further investigation, as noted above. 

4.1.1. Comparing MARCS with ISO-SWS 

Unfortunately, the first three discrepancies coincide with 
wavelength regions not free of calibration problems: the 
first one corresponds to that affected by presaturation. 
We are however confident that this issue was dealt with 
in an appropriate way as discussed in Sect. 12.11 and that 
the discrepancy seen is of a physical nature. As to the 
2nd region: this is the blue side of band 2a, which shows 
strong memory effects as shown in Sect. 12.11 and Fig. [21 
The same is true for the 3rd region, but the discrepancy is 
much wider than the wavelength band affected by memory 
effects (cf. Fig. |2J . Since a correction for memory effects 
was applied, we can at this point not rule out the possibil- 
ity that these are also true discrepancies between model 
and reality. 

From Fig. [3] and H3 it is clear that the first and fourth 
discrepancy, i.e. the strongest ones, can be due to unpre- 
dicted absorption by water, while the second one appears 
to correspond to SiO. Some contribution of extra CO ab- 
sorption might also be present in the 2.38 and 5 /xm dis- 
crepancies. 

However, for the parameters being studied, the model 
predicts almost no absorption by water: the atmospheric 
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Fig. 6. A comparison between the ISO-SWS-spectrum 
(grey) and our synthetic MARCS spectrum (black). A 
blackbody at T=3600K is shown by the dotted line for 
comparison. Clearly, our model can reproduce the global 
spectral shape very well. Nevertheless, there are some in- 
teresting discrepancies. 



temperature is too high, except for the outermost lay- 
ers. Lowering the effective temperature is not an option, 
because it is well constrain ed by the apparent diameter 
and the bolometric flux fe.g. lPerrin et al.ll2004a|) . Using a 
higher oxygen abundance does not work either since a very 
large increase would be required, which then causes far too 
strong spectral features by OH and SiO. Also an attempt 
at increasing the H2O absorption with a higher surface 
gravity fails: a physically impossible value of \ogg = 2 
would be required (see Sect. 13.11 for a discussion on the 
determination of the gravity). 

Concerning the excess emission by SiO at 4.2 /xm: since 
the SiO band at 8 /xm is well predicted by the model, it is 
most likely only a memory effect and therefore not real. 

We conclude that we see clear extra absorption by 
H2O which cannot be predicted by the hydrostatic MARCS 
model with reasonable input parameters. The problem 
thus appears not to be in the input parameters but in 
the basic assumptions. Three of these are not very likely 
to hold in the atmosphere of a variable supergiant: LTE, 
homogeneity and hydrostatic equilibrium. 

LTE is most definitely violated in the outermost 
layers of the photosphere where the density is low 
(see e.g. iRvde et al.l l2f)02(l . However, non-LTE effects 
are generally only seen in high-resolution spectra (e.g. 
lAvres fc WiedemaniJl 989L and are not believed to create 
pseudo-continuous effects such as seen here in a Orionis. 

The second assumption, homogeneity, is probably also 
violat ed in the atmosphere of Betelgeuse. iFrevtae et alJ 
(2002) calculated convection models for Betelgeuse and 
found large convective cells with strong temperature gra- 
dients between the border regions and the center. Possibly, 
more molecules, including H2O, could form there too. 
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The last assumption has been well studied in the con- 
text of Mira stars. Hydrodynamic models are essential in 
the modelling of observed IR spectra of these pulsators 
fWoitke et al.ll999llHofner et all200.4 . Although the am- 
plitude of the visual variations of a Orionis is small com- 
pa red to that of a typical Mira, i.e. abo ut 0.25 mag in 
V dGravl^OOof) vs. 3 mag for e.g. T Cep l|Van Malderenl 
120031) . the low surface gravity and high luminosity may 
help the levitation of the upper layers. As for Mira stars, 
the temperatures in this levitated matter can be low 
enough for molecules to form. If shocks are formed, by 
outward moving material running into infalling material, 
high molecular density shells will result. 

Extra molecular layers are thus an attractive solution 
for the observed discrepancies. We will investigate this 
hypothesis below, but let us first remark that the an- 
gular diameter on the sky required to fit the ISO-SWS 
spectrum is only 43.6 mas (limb-darkened diameter of the 
r ross=10 -7 layer), to be compared to an observed LD di- 
ameter of 45.6 mas (obtained by fitting L band visibili- 
ties computed from our MARCS model to the TISIS ob- „ 
servations 5 ). This confirms the hypothesis formulated in ^ 
Sect. EIH that the absolute calibration of the ISO-SWS 
spectrum based on band Id has caused an underestimation 
of the actual flux, resulting in both the large deviations 
from unity for the scaling factors presented in Tableland 
the discrepancy in photometric and interferometric diam- 
eter reported here. The interferometric diameter would 
suggest the actual flux to be about 10% higher than pre- 
dicted by band Id (cfr. Tabled). 

4.1.2. Adding an H 2 layer 

Fitting the ISO-SWS spectrum with a 
photosphere+water-layer model, by a \ 2 minimiza- 
tion in the grid presented in Table |31 we find the best fit 
for a layer of 1750 K, with an outer radius of 1.45 R* and 
an H2O column density of 2 x 10 19 cm~ 2 . 

A few discrepancies remain, and they can be identified 
as due to CO (at 2.4 and 4.5 /xm) and OH (several strong- 
lines between 3 and 3.5 /xm). The fit to the ISO-SWS spec- 
trum can be improved some more by adding 2 x 10 21 cm~ 2 
of CO and 5 x 10 20 cm~ 2 of OH. The resulting synthetic 
spectrum is shown in Fig. [7] 

4.1.3. Comparison with the near-IR interferometric 
data 

From this model, we can compute K and L wide band visi- 
bilities, to be compared to the FLUOR (K) and TISIS (L) 
observations. The result is shown in Fig. IH1 The agreement 
is convincing for a photospheric limb-darkened diameter 
of 45.6 mas. In fact, the layer opacity in the FLUOR and 
TISIS bandpasses is so low that there is no noticeable dif- 
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Fig. 7. The ISO-SWS spectrum (black) is compared to 
the MARCS photosphere model (green), and the MARCS 
photosphere model embedded in a molecular layer (red) 
at 1.45 R*, with a temperature of 1750 K and H2O, CO 
and OH column densities of 2 x 10 19 cm~ 2 , 2 x 10 21 cm~ 2 , 
and 5 x 10 20 cm~ 2 respectively. The improvement is signif- 
icant. The problems at 4.2 and 5 /zm can be attributed to 
memory effects in the ISO-SWS observation (see Fig. 01. 



5 The agreement between photosphere model and ISO-SWS 
spectrum is excellent across the L band. There should be no 
contribution of extra layers at these wavelengths 



ference in visibility curve between pure photosphere and 
photosphere+layer model. The derived photospheric di- 
ameter i s slightly lar g er than determined from the same 
data bv IPerrin et alJ l)2004al) . but this is due to the ex- 
tent of our stellar photosphere model: the outermost layer 
of our model corresponds to R rrosB= io-7. This translates 
into a "tail" on the intensity pr ofile which is not pr esent 
in the analytical profile used bv lPerrin et alJ 1 2004af) and 
which amounts to about 4 percent of the total stellar di- 
ameter. We conclude that the photospheric diameter of 
our best-fitting photosphere +layer model is in go od agree- 
ment with the one found bv lPerrin et alJ l|2004a|) . i.e. that 
the molecular layer does not change the apparent size in 
the near-IR at low spectral resolution. 
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Table 3. The grid of layer parameters in which we search 
for a best fit to both the ISO-SWS spectrum and the in- 
terferometric data. 



parameter 




min. val. 


max. val. 


stepsize 


Rout [R*] 




1.20 


1.50 


0.05 


Temp. [K] 




1500 


2500 


100 


Col. dens. H 2 


[cm" 2 ] 


6 x 10 18 


4 x 10 20 


x 2 


Col. dens. CO 


cm" 2 ] 


1 x 10 20 


1 x 10 24 


x 10 


Col. dens. OH 


cm" 2 ] 


1 x 10 20 


1 x 10 24 


x 10 


Col. dens. SiO 


[cm- 2 ] 


1 x 10 20 


1 x 10 24 


x 10 




20 40 60 

Spatial Frequency (arcsec ) 



Fig. 8. Observed K and L band visibilities are compared 
to our best single layer model. The match is excellent for 
a photosphcric LD diameter of 45.6 mas. The very dis- 
crepant points at low spatial frequency are TISIS L band 
observations, which might be corrupted by a poor sub- 
traction of the thermal background. 



4.2. The mid IR 

4.2.1. Excess emission and a larger apparent diameter 

The mid-IR part of the ISO-SWS spectrum is dominated 
by the Si-0 stretching and O-Si-0 bending resonances 
in amorphous olivines at 9.7 and 18 /im respectively 
(Fig. [TDjl . The MIDI N-band spectra on the other hand 
do not show such a 9.7 /im feature (Fig.0. The slit used 
for the latter observations is only 0.52 arcsec wide. This 
confirms the la rge inner radius of the olivine dust shell 
as suggested bv lSloan et alJ (^)93) who found no silicate 
emission inside a region of about 0.5 arcsec, some emis- 
sion between 0.5 and 1 arcsec, and the better part of the 
silicate emission even further out. This is also consistent 
with U KIRT mid-IR images presented by [Skinner et alJ 
<ll997ft and w i th the inner radius (1 arcsec) measured by 
iDanchi et alJ l|l994f) with the ISI interferometer. 

Surprisingly, there is excess (non-silicate) emission 
within the PSF of the MIDI N-band spectra which in- 
creases with wavelength up to about 11 /im and then levels 
out (Fig. EjJ • That such excess emission is present is also 
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Fig. 9. Average MIDI N-band spectrum compared to the 
MARCS prediction for the photospheric emission. Upper 
panel: the calibrated MIDI spectrum and the MARCS pho- 
tospheric spectrum. There is no trace of olivine emission 
in the MIDI slit, which is 0.54 arcsec wide, confirming that 
the olivine dust is indeed more than 10R* out. Lower 
panel: the MIDI spectrum divided by the MARCS spec- 
trum. There is excess emission within the PSF which in- 
creases with wavelength up to about 11 — 12 /im and then 
levels out. 

Table 4. Dust shell parameters for the best fit to the 
ISO-SWS spectrum. The dust shell inner radius corre- 
spond s to the measured inne r radius of ISkinner et alJ 
l|l997|) and lSloan et alJ lll993l) The dust to gas ratio is 
taken from Knapp et al. and the outflow velocity 

from lKnapp fc Morris! 1)19851) . Grain shapes are a contin- 
uous distribution of ellipsoids (CDE). 



Parameter 



Value 



Mdust 
Voo 

Dust to gas ratio 

Rin 
Rout 

Composition 
Grain size 
Grain shapes 



6.3 x 10- 10 M Q yr- 1 
15 km s" 1 
0.0025 
20 R* 

oo 

Amorphous MgFeSi04 
0.01 - 1 /im 
CDE 



confirmed by a MODUST modelling of the dust shell based 
on the ISO-SWS spectrum: in Fig. ^3 we see that the 2 
features due to amorphous olivine can be modelled quite 
accurately but that the flux predicted in between is far too 
low. This problem cannot be solved with other grain sizes, 
nor with another extent of the dust shell. However, when 
subtracting the excess emission seen in the MIDI spec- 
trum from the ISO-SWS spectrum, we arrive at a far bet- 
ter agreement between model and observations. Moreover, 
the excess seems to decrease again beyond 13 /im, disap- 
pearing entirely at 17 /im. 
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Fig. 10. The ISO-SWS spectrum of a Orionis (black) in 
the mid-IR is dominated by the amorphous olivine reso- 
nances at 9.7 and 18 /im respectively. If the dust features 
are modeled with an olivine outflow (parameters listed in 
Table |3J, a large discrepancy from 8 — 17 /im suggests 
another source of emission/opacity. Indeed, when we sub- 
tract from the ISO-SWS spectrum the excess emission 
found in the MIDI N-band spectrum, the agreement is 
far better (solid line). 



We remark that the dust mass loss rate we de- 
rive from a best fit to the ISO-SWS spectrum (6.3 x 
10 ~ 10 yr -1 ) is in good agreement w i th th a t found 
bvlKnapp et alJ lll980UKnapp fe Morris! l|l985|) . iKnappl 
( 1986) and lSkinner k. Whitmord (1987). It is worth stress- 
ing that in the case of this optically thin dust shell, it is 
important to take into account the photospheric SiO band 
head at 8 /im when using the 9.7 /jm feature to determine 
the dust mass loss rate. 



IWeiner et all (2£ 
11.15 /im) performed 



present ISI observations (at 
between November 1999 and 
December 2001 at spatial frequencies which cover very 
well the first lobe of the visibility curve due to the central 
object, i.e. the object inside e.g. the MIDI PSF. Together 
with earlier observations at much smaller baselines which 
allow the determination of flux ratio between central ob- 
ject and olivine dust shell, these allow a diameter deter- 
mination of the central object at that wavelength. The 
flux ratio determi ned from these observations ranges from 
55 - 65 percent l|Danchi et all Il994t iBester et all 1199(4 
ISudol et al .I F 9991 which is compatible with the flux from 
the central object in the MIDI spectra. The diameter is 
about 55 — 60 mas l|Weiner et al.ll2003|) . i.e. almost 1.5 
times the size in the near IR. Moreover, no data in the 
second lobe of the visibility curve are available, and limb- 
darkening may therefore be significant, increasing the size 
of the central object at 11.15 /im even more. 



4.2.2. Explaining the mid-IR data 

In the previous section, we found that in the mid IR, 
a Orionis (without the detached olivine dust shell) ap- 
pears to be about 1.5 times as large as in the near IR. 
Moreover, we found excess emission from 8.5 /im onward, 
which increases smoothly with wavelength up to about 11 
or 12 /im, where the excess levels out at 25 — 30 percent. 
Beyond 17 /im there appears to be no longer an excess. 
There are 3 possible sources of extra opacity which could 
be responsible for the observed excess and diameter in- 
crease: (1) an extra molecular layer, (2) chromospheric 
emission and (3) dust. Below, we investigate each of these 
possibilities. 



An extra molecular layer From an analysis of the near-IR 
part of the ISO-SWS spectrum, we found evidence for an 
extra molecular layer at about 0.45 R* above the pho- 
tosphere, containing mainly water. This putative layer is 
so optically thin that it does not influence the apparent 
size in the K and L bands. S ince H?Q opacity is in fact 
even lower around 11 /im fe.g. lVan Malderenl2003h . we do 
not expect a significant diameter increase at those wave- 
lengths. Indeed, to reproduce the ISI observations, the 
layer has to contain about 2 x 10 22 cm -2 of water (up- 
per panel of Fig. llip . which is 1000 times more than what 
we found from the near-IR part of the spectrum. Such a 
thick H2O layer is clearly not present: it would create op- 
tically thick water bands in the near-IR as well, and those 
are not seen in the ISO-SWS spectrum. Moreover, H 2 
opacity has a minimum around 10 /im which would yield 
a clear spectral signature in the MIDI N-band spectrum 
(unless the layer is extremely thick and isothermal), which 
is not observed. We conclude that H2O opacity can not be 
responsible for the observed diameter increase and excess 
emission. 

However, there may be other molecular sources of 
opacity at 11.15 /im. As can be seen from Fig. [31 both SiO 
and OH show spectral lines around 11 /im (TiO, which is 
not shown in Fig. does not have significant lines in the 
mid-IR) . Those of OH are well separated an d not very nu- 
merous: our line list ((Goldman et al .ill 998^ has no strong 
lines in the ISI bandpass. The red wing of the opacity pro- 
file of hot SiO might cover the ISI bandpass, but a very 
large column density would be required to reach the re- 
quired amount of opacity. This is not compatible with the 
lack of an extra-photospheric SiO signature around 8 /im 
in the ISO-SWS spectrum. The lower panel of Fig. ^dis- 
plays the opacity of H2O and SiO in the ISI bandpass. We 
thus exclude a molecular cause for the diameter increase 
and excess emission is. 



Chromospheric emission From UV and radio observations, 
a Orionis is known to have both a hot chromospheric com- 
ponent and an ionized wind with fairly low electron tem- 
peratures (e.g .|Basri et al.ll98lilGilliland fc Dupreefl l996: 
il.l ll997ULim et alJll998h . lHarper et all l(200lh 
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Fig. 11. Upper panel: Comparison of different models 
with the ISI obs e rvatio ns (diamonds) at 11.15 /im made 
by IWeiner et alJ l|2003|b To match these data with an 
H2O layer, we need a high H2O column density of about 
2 x 10 22 cm -2 , which is neither compatible with the near- 
IR spectrum, nor with the MIDI N-band spectrum. Lower 
panel: Opacity in the ISI bandpass (shaded region) of 
2 x 10 19 cm" 2 of H 2 (solid line) and 2 x 10 21 cm" 2 of SiO 
(dotted line), both at a temperature of 1750 K. The opac- 
ity profile is redshifte d by 21 kms -1 , which is the radial 
velocity of a Orionis l|Wilsonlll95i . 



present a model of the gas in the extended atmosphere of 
a Orionis from 1 — 10 R*, which matches the radio data 
up to the mid IR. It consists of a weakly ionized wind with 
mainly H _ f-f opacity. The temperature distribution peaks 
at about 80 mas (diameter) or 1.45 times the (mid-IR) 
ISI diameter. The maximum mean electron temperature 
in their model is only about 4000 K, but since that is not 
enough to explain the UV observations, they suggest an in- 
homogeneous medium: the region around 1.45 Rigi has an 
ambient temperature of only 2000 K, but contains a mul- 
titude of chromospheric hot spots with electron tempera- 
tures of the order of 10000 K. With an areal filling factor 
(AF) of less than 0.25, this could reconcile the UV traces 
of a hot chromospheric component with a mean electron 



temperature of only 4000 K. While the mean model pre- 
dicts neither excess emission nor a larger apparent size at 
11.15 /im (Harper, private communication), we investigate 
whether the suggested chromospheric hot spots could be 
at the origin of the mid-IR excess and diameter increase. 

We have strong constraints on the spectral shape of 
the excess emission and opacity: it should be transparent 
up to 8.5 /j,m, then increase strongly up to about 12 /im 
and decrease again toward longer wavelengths (Fig. Inland 
I10L In the optically thin regime, the excess due to chro- 
mospheric emission is constant throughout the N band be- 
cause the free-free opacity is proportional to A 2 and the 
chromospheric temperature is high enough for the source 
function to be in the Rayleigh-Jeans limit. Such a con- 
stant excess is not compatible with the observed increase 
of the excess with wavelength (see also Fig. II 3f) . In the 
optically thick regime, the observed shape of the excess 
requires (1) the temperature of the chromosphere to de- 
crease strongly along the line of sight and (2) the opacity 
to reach unity within the N band. Because the total ex- 
cess is far below that of a homogeneous, hot and optically 
thick chromosphere with a diameter of about 1.5 R*, the 
chromosphere must be very clumpy. 

To test this hypothesis, we constructed a model mim- 
icking the central photospheric disk embedded in a halo 
with optically-thick chromospheric hot spots. This model 
yields the contrast ratio between the mean 6 observed pho- 
tospheric intensity and mean layer intensity as a function 
of layer diameter and chromospheric AF factor, under the 
constraint that the chromospheric hot spots must generate 
the excess emission seen at 11 /an in the MIDI spectrum 
and the visibilities observed with the ISI. Fig. ^] shows 
how each assumed diameter of the chromospheric layer re- 
quires a certain AF factor and hot spot temperature, if it 
is to match the observed excess and visibilities. The AF 
factor goes down asymptotically to zero for a layer radius 
approaching 1.8 R*. At the same time, the chromospheric 
temperature increases to infinity, but it is clear that our 
model is too crude an approximation of the actual inten- 
sity distribution on the sky in this regime. 

We find that for almost the entire range of possible 
layer diameters, the hot spot temperature is far below the 
temperature of an actual hot chromospheric component 
(10000 K). Only in the asymptotic regime toward a layer 
radius of 1.8 R* (80 mas on the sky) does the model predict 
such a high temperature, at a very low AF factor (AF 
< 0.03), but the agreement with the observed visibilities 
is not as good as with a smaller layer diameter. Moreover, 
it appears unlikely that the very low AF factor can be 
combined with the requirements on the source function. 

6 Note that the approximation of the spotted intensity profile 
with a mean intensity profile yields reasonable visibility curves 
only if the size of the spots is much smaller than the spatial 
resolving power of the observations. This assumption is not vi- 
olated because (1) the observed visibilities are all within the 
first lobe, and (2) they are all compatible with a single layer di- 
ameter and therefore argue against large scale deviations from 
spherical symmetry at these wavelengths. 
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Fig. 12. Both the visibilities and the excess emission at 
11 jLtm can be explained as due to a halo covered with 
hot spots. Allowed layer diameters range from about 1.3 
to 1.8 R*. The derived hot spot temperatures remain well 
below the expected chromospheric temperature (10000 K, 
er et afl l200l|) . 



Harper 



We conclude that, while there must be some emission 
by a hot chromospheric component in the mid IR, the 
hot spot temperatures and required areal filling factor, 
together with the spectral shape of chromospheric emis- 
sion, appear to be in contradiction with the hypothesis 
that an inhomogeneous hot chromosphere is responsible 
for the observed N-band excess and diameter increase. 



Dust The last viable candidate source for the extra opac- 
ity and emission is dust. The maximum of the excess emis- 
sion appears to be located close to 13 /mi, which suggests 
a relation with the "13 /im feature" as seen in about 50 



JSloan et al 


1996) 


IPosch et al. 


1999), 



Si02 ( Sp eck et al.1 

i . 7i ., 



2000) 



1989). In fact, all but the last candidate have a fairly nar- 
row, peaked absorption coefficient which can not explain 
our observations. 

A simultaneous modelling of the N-band excess and the 
ISI visibilities with a layer of AI2O3 results in a very good 
agreement (Fig. I13fl for a layer at 1900 K and a column 
density of 1 x 10~ 3 g cm -2 . Moreover, the transparency 
of AI2O3 in the near IR makes even the SiO band head at 
7.7 /im well visible through the layer, as required by the 
ISO-SWS spectrum. 

The derived temperature should be confronted with 
the temperature regime only 0.5 R* above the stellar pho- 
tosphere. Although the effective temperature of a Orionis 
is 3600 K, the outermost layers of our MARCS model have 
temperatures of the order of only 2000 K. It appears there- 
fore not unlikely that the region 0.5 R* above the photo- 
sphere has a gas temperature of about 1900 K. However, 



„ 6000 



MIDI spectrum 
MARCS photosphere 
MARCS + Al 2 3 at 1900K 
MARCS + Al 2 3 at 1 550K 
MARCS + chromosphere 




Wavelength [/^m] 

Fig. 13. By including 1 x 10~ 3 g cm~ 2 of amorphous 
AI2O3 at 1900 K in the layer at 1.45 R*, we find a good 
fit to the ISI data (diamonds in the upper panel), and 
to the N-band excess (lower panel). The minor but sys- 
tematic deviation at small spatial frequency suggests that 
the layer is actually located a little closer to the photo- 
sphere, i.e. 1.35 R*. The grey box indicates the flux excess 
allowed by the flux ratio between the detached dust shell 
and the central object, i.e. the scaling factor used on the 
interferometric observations in the upper panel. The ex- 
cess of an optically thin is othermal chromosphere as in 
ISkinner fc Whitmord 1^987) is also shown (dotted line). 



demanding radiative equilibrium 7 , i.e. the absorbed stel- 
lar radiation should be emitted thermally, we arrive at a 
much higher temperature, of the order of 2400 K. 

AI2O3 dust grains are believed to condense in chemi- 
cal equilibrium at a temperature of 1900 K only in high- 
pressure environments, i.e. for a total pressure above 
10- 2 bar (fodders fc Feglevll 99Qh . This required pressure 
is a factor of 10 4 larger than the pressure in the outermost 
layers of our MARCS photosphere. 



7 We used the dust radiative transfer code modust and 
the optical constants of amorphous AI2O3 measured by 
iBeeemann et alJ Jl997l) and lKoike et alJ <199Sl) 
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Fig. 14. Light curve of a Orionis at the time of the ob- 
servations presented in this paper. Constructed from ob- 
servations by the American Association of Variable Star 
Observers (AAVSO) (Waagen, E. O. 2004, private com- 
munication) 



Which one of the three? A molecular origin can be ruled 
out right-away as long as we assume that there have 
not been drastic changes in the molecular environment 
of a Orionis between the epoch of the ISO-SWS and the 
near-IR interferometric observations and that of the mid- 
IR observations. The light curve of a Orionis (Fig. I14|l 
shows no peculiar behaviour in that time span. 

Chromospheric emission is certainly present, and we 
cannot rule out some influence in the mid IR. However, as 
discussed above, we see no way of reconciling the spectral 
and spatial constraints with a chromospheric origin. 

Dust appears to be the most attractive solution: it 
can reproduce spectral and spatial observations simul- 
taneously. Nevertheless, the required temperature and 
pressure are a reason for concern. Furthermore, this in- 
evitably leads to questions on the relation between such 
an Alumina layer and the detached olivine dust shell. 

4.3. Discussion 

The water layer Let us start this discussion by returning to 
the near-IR observations. We found evidence for an opti- 
cally thin molecular layer containing H2O and maybe also 
CO and OH about 0.45 R* above the photosphere. The 
derived column density (about 2 x 10 19 cm" 2 ) is slightl y 
higher than that derived bv Ijennings &: Sadal (Il998h. 
but low er than that from iTsuiil fcOOObh and lOhnakal 
ll2004bh. as listed in Ta ble The lower value found by 
Ijennings k, Sadal l)l998|) is the consequence of a higher 
adopted water temperature: 2500 — 2800 K. However, this 
temperature is o nly an upper limit on th e actual t e mpera- 
ture of the H 2 \ Jennings k Sadall998|) . l0hnakal l|2004bj) 
on the other hand infers both a higher temperature and 
a significantly larger column density. This is the conse- 
quence of his attempt to explain both the mid-IR water 
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Fig. 15. ISO-SWS spectrum (black), MARCS photosphere 
model (green) and a ph otosphere - !- layer m odel (red) using 
the layer parameters bv lOhnakal l|2004bl) . 

spectrum at 6 and 11 fim and the ISI diameter at 11.15 fim 
as being due solely to a layer of HTO^Fig. 1151 compares a 
model using the layer parameters of Ohnaka (2004b) with 
the ISO-SWS spectrum over the full wavelengt h range. 
Improv ements w.r.t. the model as presented by lOhnakal 
(2004tJ) are the use of a MARCS photosphere instead of 
a blackbody approximation and the use of the Ames line 
list for H2O instead of that from the hitemp database. 
While this model reproduces quite well the spec t ral fea- 
tures between 6 and 7 fim also used by lOhnakal §004b) 
(see inset of Fig. 1151 . it is not compatible with the spectro- 
photometric flux levels. Furthermore, the predicted band 
strength at 2.8 fim is too large and the SiO band head at 
8 /im, clearly seen in the ISO-SWS spectrum, is no longer 
visible through such a thick layer. We conclude that the 
model as presented by lOhnakal (j2004b) has too large an 
H2O column density. 

The dust layer From the mid IR excess and the 11.15 fim 
intcrferometry, we found evidence for the presence of hot 
amorphous AI2O3 at about the same location above the 
photosphere as the H2O layer. This is at least a factor of 
10 closer to the star than the onset of silicate emission. 

Several dust condensation scenarios ascribe a crucial 
role to AI2O3: it is assumed to be the first dust species 
to condense, and thereafter act as a seed nucleus for fur- 
ther dust condensation (e.g. ISalpeterlll977l: lTielensl ll990: 
Sedlmavrl Il997t lLodders k Feglevl ll999|T Remark that 
Patzen(|2004^ argue that this scenario is only valid in 
chemical equilibrium conditions, and therefore not in the 
rapidly expanding wind, but the conditions in the sta- 
tionary layer we consider here may be more favourable to 
attain chemical equilibrium. 

If indeed we have observed the very onset of dust for- 
mation, then it is quite puzzling why no dust is seen in 
between 1.5 and 10—20 R*. If the mass loss is not episodic, 
there are only two possibilities: either (1) there is no dust 



14 



Verhoelst et al.: Amorphous alumina in the extended atmosphere of a Orionis 



Table 5. St ellar and layer param eters for a Orionis, as derived bv lJennines & Sadal (|l998^ . iTsuii l|2000rJ) . lOhnakal 
(2004b]) and IPerrin et alJ l)2004a|) . The first 2 determinations are purely spectroscopic, the third combines high- 
resolution spectroscopy and interferometry and the fourth one combines photometry with interferometry. (NA = 
not available) 



Parameter Jennings fc Sada (1998) Tsuii (2000b) Ohnaka (2004b) Pcrrin ct al. (2004a) This work 



9+ [mas] 


NA 


NA 


44.0 


42.00 ± 0.06 


45.6 


T* [K] 


NA 


3600 


3600 


3690 ± 50 


3600 


flayer [mas] 


NA 


NA 


63.8 ±6.4 


55.78 ± 0.04 


65 


Tl a y cr [K] 


< 2800 


1500 ± 500 


2050 ± 100 


2055 ± 25 


1750 


Layer Composition 


H 2 


H 2 


H 2 


undefined 


H 2 0, AlaO s (?) 


K, opacity 


NA 


NA 


NA 


0.060 ± 0.003 


NA 


L opacity 


NA 


NA 


NA 


0.026 ± 0.002 


NA 


11/^m opacity 


NA 


NA 


NA 


2.33 ±0.23 


NA 


H 2 Q col. d. [cm" 2 ] 


3±2 x 10 18 


~ 1 x 10 20 


1 - 4 x 10 20 


NA 


2 x 10 19 



in this region, in which case it must be destroyed right af- 
ter its formation, or (2) there is an outflow of amorphous 
AI2O3 but at such a low density and temperature that we 
do not see it. 

The former scenario is plausible given the presence of 
a patchy hot chromospheric component which could de- 
stroy the Alumina before it reaches cooler regions. The 
dust seen at larger distances may well be formed anew, 
if the outflowing gas becomes again cool enough for con- 
densation to occur. The silicate emission then does come 
from a detached shell. 

Th e second hypothesis is supported by lOnaka et alJ 
l)l989t) . who propose that the alumina grains are highly 
transparent up to the point where they collect silicates on 
their surface. Coupled with a very low density due to a 
strong acceleration at their birth in the molecular layer, 
this would make them invisible up to the silicate condensa- 
tion location, where silicates settle onto the AI2O3 grains. 
We computed radiative equilibrium temperatures for both 
silicates and AI2O3 at the radius of the silicate condensa- 
tion and find the Alumina to be cooler by a few hundred 
K, making it indeed undetectable 8 ! 

This could mean that the dust shell of a Orionis is 
not really "detached" as previously assumed, but rather a 
continuous outflow from close to the stellar photosphere, 
which is transparent up to the point where silicates con- 
dense onto the alumina grains. Pure AI2O3 is then only 
visible at fairly high temperatures near the stellar surface 
where it is formed. 

However, a major shortcoming in this scenario is that 
Alumina is fairly transparent at short wavelengths and 
therefore radiation pressure on a 0.01 /im grain is insuffi- 
cient (by a factor of 10 at least) to initiate the outflow. 



8 For the Alumina at the silicate condensation radius to 
be detectable in the ISO-SWS spectrum, we would need a 
Alumina mass loss rate which is 10 times higher than the sili- 
cate mass loss. 



5. Conclusions and prospects 

We have modelled the molecular (and possibly dusty) close 
environment of the late-type supergiant a Orionis. We 
took into account both spectral and spatial information 
from the near to mid-IR. The improvements over previous 
modelling attempts for a Orionis are the use of a sophis- 
ticated MARCS model for the central star, the computa- 
tion of radiative transfer through the molecular layers in 
spherical geometry, up-to-date line lists, and a larger set 
of observational constraints. 

We find evidence for an optically thin layer of wa- 
ter close above the photosphere. This layer gives rise to 
some spectral signature, but does not increase the ap- 
parent size in the near-IR w.r.t. that of the pure pho- 
tosphere. However, in the mid-IR, we find excess emission 
by amorphous silicates far out in the stellar wind (at least 
20 R* from the stellar surface) and another source of ex- 
cess emission much closer to the photosphere. The extra 
source of opacity close to the star is so optically thick 
that it increases the apparent size of the star with a fac- 
tor of 1.5 from the near- to the mid-IR. It must however 
be fully transparent up to 8 (im, and the excess emission 
appears to decrease again beyond 15 fim. We show that 
it cannot be of a molecular origin since that would induce 
strong spectral features in the near IR. Chromospheric 
opacity/emission is most definitely present at radio and 
UV wavelengths, but we see no way to reconcile the spec- 
tral and spatial properties (inhomogeneity) of a chromo- 
sphere with the near and mid-IR observations. Dust grains 
of amorphous alumina (AI2O3) do yield a good spectral 
and spatial fit, but at an uncomfortably high temperature 
(1900 K). Nevertheless, this hypothesis fits in recent dust 
condensation scenarios and we believe it to be the most 
likely solution. 

New MIDI observations at 5 different baselines, high 
spectral resolution and with simultaneous photometry are 
planned for autumn 2005 and, together with requested 
VISIR high spectral and spatial resolution observations in 
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the Q-band, will undoubtedly help to select among the 
hypotheses. 
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